Superconductivity in the 2D Hubbard model: Electron doping is different 
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A variational Monte Carlo calculation is used for studying the ground state of the two-dimensional 
Hubbard model, including hopping between both nearest and next-nearest neighbor sites. Super- 
conductivity with d-wave symmetry is found to be restricted to densities where the Fermi surface 
crosses the magnetic zone boundary. The condensate energy is much larger for hole doping than for 
electron doping. Superconductivity is kinetic energy driven for hole doping, but potential energy 
driven for electron doping. Our findings agree surprisingly well with experimental data for layered 
cuprates, both for electron- and hole-doped materials. 

PACS numbers: 71.10.Fd,74.20.Mn,74.72.-h 



It is widely accepted that the layered cuprates are 
doped Mott insulators, which are well described by the 
two-dimensional Hubbard model, at least regarding the 
insulating antiferromagnetic state of the parent com- 
pounds and the transition to a metallic state upon dop- 
ing. Whether this model embodies superconductivity 
has been debated since two decades, but recent results 
obtained with improved numerical techniques strengthen 
the case of pairing induced by on-site repulsion. Progress 
has been made both in dynamical mean-field theory 
where an effective attraction is deduced from the irre- 
ducible two-particle vertex, and in variational calcula- 
tions [3| , where a broken-symmetry ground state with 
a d-wave superconducting order parameter is found in a 
certain density range and for large enough values of U. 

The Hubbard Hamiltonian H ~ Hq + UD consists of 
a hopping term ( "kinetic energy" ) 



Ho 



(1) 



and an on-site repulsion UD, where D 
the number of doubly occupied sites, n^, 



and 



the operator c^^ (cio-) creates (annihilates) an electron at 
site i with spin a. 

In our previous work [J we have restricted ourselves 
to nearest-neighbor hopping (i^ = t for nearest-neighbor 
sites, otherwise). We have used the variational ansatz 



(2) 



which yields a substantial improvement with respect to 
the Gutzwiller wave function {h = 0), and comes very 
close to the exact ground state both for small 2D sys- 
tems and for the solvable 1/r chain Both antifer- 
romagnetic and superconducting states have been chosen 
as mean- field states |^'o)- Antiferromagnetism was found 
to prevail at half filling, while d-wave pairing has been 
obtained for the doped system, below a hole concentra- 
tion of 0.18. At first sight these results seem to contradict 
recent work by Aimi and Imada Q , who see no signature 
for superconductivity in their Monte Carlo simulation. 
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FIG. 1: (Color online) Fermi surface for three different elec- 
tron densities. 



A closer examination shows that there is no discrepancy. 
On the one hand, the hole densities considered by Aimi 
and Imada for up to 10x10 lattices (0.18, 0.22, 0.28) are 
in a region where we obtained a vanishing superconduct- 
ing order parameter On the other hand, for the only 
other density studied by Aimi and Imada (0.09) their re- 
sult is inconclusive. In fact, for their parameter values 
(U/t=6, lattice size 8x8) the order parameter is expected 
to be too small to lead to a contradiction 0]. 

With hopping restricted to nearest-neighbor sites the 
Hubbard Hamiltonian (on the square lattice) is electron- 
hole symmetric and there is no difference between elec- 
tron and hole doping. However, this restriction leads to 
a (bare) Fermi surface which disagrees qualitatively with 
photoemission experiments on layered cuprates (8[. In 
the present paper we investigate the more realistic case 
where hopping between both nearest {t) and next-nearest 
neighbors {t') is included. We use parameters U = 8t and 
t' = —0.3t throughout. The bare single-particle spectrum 

e^: = — 2t(cos kx + cos ky) — At' cos k^ cos ky 

leads to the Fermi surfaces of Fig. [1] The innermost line 
corresponds to the van Hove filling where the Fermi sur- 
face passes through the saddle points at (tt, 0) and (0, tt). 
These crossings between the Fermi surface and the mag- 
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n 




g 


h 


E/t 


0.7500 


-0.9921(1) 


4.2(1) 


0.113(2) 


-0.858(1) 


0.7800 


-0.9612(3) 


4.0(1) 


0.112(2) 


-0.829(1) 


0.8125 


-0.9107(3) 


3.9(1) 


0.111(2) 


-0.795(1) 


0.8400 


-0.788(1) 


3.7(1) 


0.111(2) 


-0.763(1) 


0.9000 


-0.728(1) 


3.8(1) 


0.111(2) 


-0.676(1) 


0.9500 


-0.603(1) 


4.0(1) 


0.114(2) 


-0.591(1) 



TABLE I: Chemical potential, parameters g and h and energy 
per site for hole doping (n < 1) and an 8 x 8 lattice. 



n 




g 


h 


E/t 


1.0500 


0.7666(1) 


3.6(1) 


0.109(2) 


-0.222(1) 


1.0800 


0.6440(1) 


3.4(1) 


0.106(2) 


-0.069(1) 


1.1000 


0.5488(1) 


3.2(1) 


0.104(2) 


0.040(1) 


1.1300 


0.4380(1) 


3.0(1) 


0.100(2) 


0.206(1) 


1.1600 


0.3870(1) 


2.9(2) 


0.096(3) 


0.374(1) 


1.2000 


-0.2996(1) 


2.6(2) 


0.091(3) 


0.608(1) 



TABLE IL Chemical potential, parameters g and h and en- 
ergy per site for electron doping (n > 1) and an 8 x 8 lattice. 



netic zone boundary (the "hot spots") move inwards as 
the density n is increased and finally merge (outermost 
line). Hot spots arc restricted to 0.726 < n < 1.206. Our 
results, to be discussed below, indicate that supercon- 
ductivity occurs only in this range. 

Here we use again the variational ansatz of Eq. ([T]) , but 
restrict ourselves on a d-wave superconducting ground 
state, which introduces, in addition to g and h, two 
other parameters, the amplitude A of the superconduct- 
ing gap function and the "chemical potential" ^. To com- 
pute the variational energy, the exponent of the operator 
e^f^ is first decoupled by applying a discrete Hubbard- 
Stratonovich transformation [9|, which introduces an 
Ising spin at each site. All operators are then quadratic 
in creation and annihilation operators and therefore the 
fermionic degrees of freedom can be integrated out. The 
remaining sum over Ising spin configurations is performed 
by a Monte Carlo simulation. In order to avoid the minus 
sign problem, calculations are carried out in the grand 
canonical ensemble with an average density fixed by /i, 
which is therefore not a variational parameter. To reduce 
the statistical error, the optimization procedure is based 

13, El. 



We 



on the method proposed by Ceperley et al. 
have used periodic-antiperiodic boundary conditions. 

The minimization of the energy E for fixed average 
densities yields the results of Tables I and II for hole 
and electron doping, respectively. We notice that the pa- 
rameter which controls the crossover between itinerant 
(small g) and "localized" (large g) many-particle states, 
remains large for hole doping, but decreases rapidly for 
electron doping. Therefore, while the hole-doped region 
0.75 < n < 0.95 is a "localized" doped Mott insulator, 
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FIG. 2: (Color online) Optimized gap parameter as a function 
of doping for two different lattice sizes. For comparison, the 
result for t' = is also shown (from Ref. [j). 




FIG. 3: (Color online) Superconducting order parameter as a 
function of doping for an 8 x 8 lattice. 



the electron-doped part 1.05 < n < 1.2 rapidly undergoes 
a crossover to an itinerant regime as n increases. We at- 
tribute this difference to the bare single-particle density 
of states at the Fermi energy, which is much larger for 
hole than for electron doping. 

Figs. [Hand [3] show, respectively, the gap parameter A 
and the superconducting order parameter $ = KCi|C^.|)| 
as functions of doping concentration x = 1 —n. The cor- 
responding results for = Q are completely electron- 
hole symmetric and reproduced only in the right pan- 
els. On the hole-doped side superconductivity exists for 
< X < 0.25 (A remains finite, but <& vanishes for 
X — > 0), i.e. in a larger region than for t' = 0. In 
contrast, on the electron-doped side the superconduct- 
ing region is reduced to —0.2 < x < —0.05. Thus we 
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electron doping the decrease in potential energy gives 
the main contribution to the condensation energy. The 
same behavior has been consistently obtained for other 
densities. The quantity 



(3) 



- 0.01 also plotted in Fig. [51 is ~ up to a minus sign - propor- 
tional to the integrated optical conductivity originating 
from intraband transitions [iJl . For the simple Hubbard 



FIG. 4: (Color online) Condensation energy per site for an 
8x8 lattice. 
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FIG. 5: (Color online) Changes in kinetic, potential and total 
energies as well as in the quantity W as functions of the gap 
parameter, for an 8 x 8 lattice. {—W is proportional to the 
oscillator strength for intraband transitions.) The relative 
uncertainties are smaller than the symbol sizes. 



find indeed that superconductivity is restricted to den- 
sities where the (bare) Fermi surface passes through hot 
spots (see Fig. [1]). The qualitative difference between 
electron- and hole-doping for x ^ suggests that super- 
conductivity competes strongly with antifcrromagnctic 
order (which we did not take into account in this study) 
on the hole-doped side and weakly on the electron-doped 
side. Further support for this conclusion comes from the 
condensation energy Econd = E{0) — E{A), which is one 
order of magnitude smaller for electron doping than for 
hole doping, as depicted in Fig. [H 

Fig. El shows the kinetic, potential and total energies 
as functions of the gap parameter for the two densities 
n = 0.78 and 1.16. For hole doping the energy gain 
is clearly due to a decrease in kinetic energy, while for 



model {t' = 0) W is equal to the kinetic energy, but for 
t' ^ the two quantities differ. For hole doping W has 
a pronounced minimum in the region of the optimal gap, 
corresponding to an increase of oscillator strength, while 
for electron doping W increases monotonically with A, 
akin to BCS behavior where the oscillator strength is re- 
duced at the onset of superconductivity. 

We attribute this asymmetry between hole and elec- 
tron doping to the different values of the correlation pa- 
rameter g (see TablesHland HI)) . which puts the hole-doped 
system into the "localized" regime, while the electron- 
doped system is more itinerant. To make the point clear 
we consider the simple Hubbard model {t' = 0) both in 
the small U (itinerant) and large U ("localized") lim- 
its. In the small U limit superconductivity is produced 
by the Kohn-Luttinger mechanism (3], where the con- 
densation energy is expected to increase as a function of 
U/t, whereas in the large U limit the condensation en- 
ergy arises from magnetic exchange and thus is likely to 
increase with t/U. The change in kinetic energy is then 
obtained through the Hellman-Feynman theorem. 



d 

Ekin{A) — Ekin{0) = —t-^Econd ■ 

at 



(4) 



which is positive in the small U limit and negative in the 
large U limit. 

Additional information can be gained from the mag- 
netic structure factor 



S{q) 



1 

TV 



)((n,| -7i,x)(njt -riji)) , (5) 



displayed in Fig. [SI S{q) exhibits a clear maximum at 
(tt, tt), which is largest close to half filling and decreases 
as doping increases. There is very little difference be- 
tween electron and hole doping, presumably because in 
the large U limit spin correlations depend mostly through 
the exchange constants J = At^ /U and J' = 4t'^/C/ on 
the microscopic parameters and therefore are essentially 
independent on the sign of t' . 

We have seen above that superconductivity is re- 
stricted to the region where two points of the Fermi 
surface can be connected (at least approximately) by 
the antifcrromagnctic wave vector (tt, tt). This together 
with the strong peak of S{q) for q = (tt, tt) supports 
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FIG. 6: (Color online) Magnetic structure factor as a function 
of the momentum, for various densities above and below half 
filling, for an 8 x 8 lattice. 



the point of view that superconductivity in the two- 
dimensional Hubbard model is due to a magnetic mech- 
anism. This conclusion is now widely accepted, but the 
question whether the mechanism is of the RVB ty pe 141 
or arises from the exch ange of spin fluctuations [l5| is 



presently under debate [la, [17|. We cannot solve this 



problem on the basis of our variational calculations, but 
the distinction between "localized" and itinerant regimes 
can give additional useful insight. According to our re- 
sults the former is appropriate for the hole-doped region, 
where superconductivity is associated with a decrease in 
kinetic energy, the latter regime is found for electron dop- 
ing, where superconductivity arises from a more conven- 
tional gain in potential energy. 

The results described above compare surprisingly well 
with experiments on cuprates, better than our earlier cal- 
culations (3] where only nearest neighbor hopping has 
been taken into account. This concerns the phase dia- 
gram, in particular for electron-doped materials for which 
a recent systematic study [l3| gives a dome-shape (for Tc) 
strikingly similar to our Fig. [2] Photoemission data give 
values of the superconducting gap A in the range 10-20 
meV for hole-doped compounds (LSCO ^ or Bi2201 
0) and ~ 5 meV for an electron-doped compound (Nd- 
CeCO [i^). Choosing t — 300 meV (neutron data), our 
maximum gap parameters are 15 and 7 meV for hole and 
electron doping, respectively. An increase of oscillator 
strength, predicted by our calculations for hole doping, 
has been reported on the basis of optical spectroscopy for 
an underdoped sample [2]| , but the situation is less clear 
on the overdoped side. We are not aware of correspond- 
ing measurements on electron-doped materials. 

In conclusion, our variational search for c?-wave su- 
perconductivity in the two-dimensional Hubbard model 
gives a differentiated picture, namely a large (moderate) 
correlation parameter for hole (electron) doping, a gain 



in kinetic (potential) energy due to pairing and a large 
(small) condensation energy. Our wave function is ex- 
pected to be better suited for describing the electron- 
doped region, which is less "localized" than the hole- 
doped region. More work is needed to improve our under- 
standing of the pseudogap phase observed in the cuprates 
for weak hole doping, by studying possible competing in- 
stabilities. 
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